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Abstract 

Oi&/e£"f/Ve5; Ventricular septal defect (VSD), one of the most common types of congenital heart disease (CHD), results from a 
combination of environmental and genetic factors. Recent studies demonstrated that microRNAs (miRNAs) are involved in 
development of CHD. This study was to characterize the expression of miRNAs that might be involved in the development 
or reflect the consequences of VSD. 

Methods: IVliRNA microarray analysis and reverse transcription-polymerase chain reaction (RT-PCR) were employed to 
determine the miRNA expression profile from 3 patients with VSD and 3 VSD-free controls. 3 target gene databases were 
employed to predict the target genes of differentially expressed miRNAs. miRNAs that were generally consensus across the 
three databases were selected and then independently validated using real time PCR in plasma samples from 20 VSD 
patients and 1 5 VSD-free controls. Target genes of validated 8 miRNAs were predicted using bioinformatic methods. 

Results: 36 differentially expressed miRNAs were found in the patients with VSD and the VSD-free controls. Compared with 
VSD-free controls, expression of 15 miRNAs were up-regulated and 21 miRNAs were downregulated in the VSD group. 15 
miRNAs were selected based on database analysis results and expression levels of 8 miRNAs were validated. The results of 
the real time PCR were consistent with those of the microarray analysis. Gene ontology analysis indicated that the top target 
genes were mainly related to cardiac right ventricle morphogenesis. NOTCHl, HANDl, ZFPIV12, and GATA3 were predicted 
as targets of hsa-let-7e-5p, hsa-miR-222-3p and hsa-miR-433. 

Conclusion:\Ne report for the first time the circulating miRNA profile for patients with VSD and showed that 7 miRNAs were 
downregulated and 1 upregulated when matched to VSD-free controls. Analysis revealed target genes involved in cardiac 
development were probably regulated by these miRNAs. 
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Introduction 

Congenital heart disease (CHD) is cardiovascular malformation 
caused by abnormal development of the fetal heart and the great 
vessels. Its incidence is 32.74/10000 in China [1] but 5.4-16.1/ 
1000 in the other countries [2,3]. CHD is one of the most common 
types of congenital malformation in children, and is also one of the 
main causes of neonatal and infant mortality [4] . 

The causes of CHD are complex, and are not fully understood. 
Most CHD (90%) may be related to genetic and environmental 
factors [5] . The common CHDs include atrial septal defect (ASD), 
ventricular septal defect (VSD), patent ductus arteriosus (PDA), 
tetralogy of Fallot (TOF), transposition of the great arteries (TGA), 
pulmonary valve atresia (PA), coarctation of the aorta (COA) and 
tricuspid atresia (TA). Among them, VSD is one of the most 
common, accounting for about 20% of CHD [6]. VSD is mainly 
caused by left and right ventricular septal defect-induced abnormal 
traffic. Despite the fact that its embryology and physiology have 
been elucidated, its etiology and pathogenesis are unclear [7]. 



A microRNA (miRNA) is a post-transcriptional regulatory 
factor and small single-stranded non-coding RNA molecule, 18— 
22 nucleotides in length. It can pair with 3' non-coding regions of 
a target gene's mRNA, and negatively regulate expression of target 
genes at the posttranscriptional level. It can regulate cell growth, 
metabolism, diflFerentiation and apoptosis, participating in the 
growth of the living organism [8] . miRNA plays an important role 
in many physiological and pathological processes. miRNA genes 
represent only 1% of human genes [9]. However, sequence 
analysis suggests that miRNA can potentially regulate 30% of 
human genes through complex regulatory networks [10]. The 
miRNA stably expressed in body fluids plays an important role in 
cardiovascular diseases and tumor occurrence, and circulating 
miRNA can be used as a potential biomarker for disease diagnosis 
as it is very stable in serum and cannot be degraded by RNA 
degrading enzymes [11-14]. Recent studies showed that miRNA 
is involved in embryonic heart development, morphogenesis of the 
heart, and myocardial cell growth and differentiation, playing an 
important function in the occurrence and development of 
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cardiovascular disease [15]. It was reported that miRNA is 
associated with the pathogenesis of CHD [16], and it plays an 
increasingly important role in diagnosis and treatment of heart- 
related diseases [17,18]. MiR-1 and miR-13.S can control 
development of the myocardium and the skeleton [19]. Huang 
et al. [20] found that miRNAs may play a central role in 
craniofacial and cardiovascular systems, and if mutated may cause 
nerve - craniofacial - congenital heart defects. A number of target 
genes regulated by miRNA that might be related with VSD have 
been identified. For example, miRNA -21 and miRNA- 181a play 
an important role in the occurrence and development of VSD in 
mice with VSD phenotype after Dicer gene knockout [20] . A very 
recent study has also found that two miRNAs, miR-1-1 and miR- 
181c, are associated with VSD pathogenesis [21]. However, 
although there are many miRNAs related to the occurrence and 
development of VSD, it is unclear how they regulate VSD, and 
more research is needed. 

miRNA microarray analysis for gene expression proliles is one 
of the most important methods used to screen and study the 
occurrence and development of disease-associated specific miR- 
NAs. We used it here to screen differentially expressed miRNAs in 
the plasma of patients with VSD and that of VSD-free controls. 
Meanwhile, quantitative real-time fluorescent polymerase chain 
reaction (RT-PCR) was used to verily the reliability of miRNA 
microarray analysis in detecting differentially expressed miRNAs. 
In addition, we predicted downstream target genes regulated by 
difierentially expressed miRNAs using target gene prediction 
software, and further analyzed its biological function. It is hoped 
that this study wUl provide valuable information for clarifying the 
pathogenesis of VSD at the molecular level. 

Patients and Methods 

Participant selection 

The study was approved by the Ethics Committee of Taishan 
Medical College and performed with the written informed consent 
of all the patients and guardians. 

Patients with VSD and VSD-free participants were selected 
from the Children's Hospital of Jinan City (Shandong, China) 
between August 2012 and June 2013. Patients with VSD were 
enrolled from surgical cardiovascular patients with CHD (shunt 
diameter of 5-8 mm) diagnosed by ultrasound technique or by 
intraoperative observation. All the patients had VSD but without 
other deformities. Meanwhile, the VSD-free participants in the 
control group were selected from hospitalized orthopedic patients. 
They were from the same areas as the patients with VSD, but had 
no CHD genetic history. One patient suffered from a fracture of 
the clavicle while the other two suffered from subluxation of the 
radial head. 

Three patients with VSD or VSD-free participants were 
selected for miRNAs microarray analysis. 20 patients with VSD 
and 15 VSD-free participants were selected to verify the 
differentially expressed miRNAs using RT-PCR. 

Blood sample collection and plasma preparation 

The fasting blood sampling was performed in the morning. 
Venous blood (5 mL) was collected from each participant in an 
ethylenediaminetetraacetic acid (EDTA) tube, allowed to stand for 
10 min, and then centrifuged at 1500 rpm for 10 min. The upper 
light yellow liquid was harvested, which was plasma. This was then 
aliquoted into RNase-free epoxy resin tubes and kept at — 80°C. 



Plasma RNA extraction and quality detection 

The total liNA in 250 j^l plasma was extracted following 
instructions of the TRIzol and miRNeasy mini kit (Qiagen, USA). 
The RNA was dc-tected by Thermo Scientific NanoDrop 1000 
Spectrophotometer and its purity was assessed by the ratio of 
absorbance (OD value = A260/A280) at 260 and 280 nm. A 
value between 1.8 and 2.1 was considered as acceptable. 

Plasma miRNA labeling and isolation 

The miRCURY Hy3/Hy5 Power labeling kit (Exiqon, 
Vedbaek, Denmark) was used for miRNA labeling. Each sample 
was 3 '-end-labeled with the Hy3 fluorescent label using RNA 
ligase following the instructions of the labeling kit to obtain 

fluorescent probes for microarray hybridization. 

The miRNA was isolated according to the kit instructions. 
Specifically, 2 |J.l RNA was dissolved in water, mixed with 1 Jll calf 
intestinal phosphatase (CIP) buffer and CIP. The mixture was 
incubated at 37°C for 30 min, and terminated by incubation at 
95°C for 5 min. Then, 3 [i\ labeling buffer, and 1.5 |tl fluorescent 
label (Hy3), 2 |a.l DMSO, and 2 ^1 labeling enzyme were added to 
the mixture. The labeling reaction was incubated at 16°C for 1 h 
and then terminated at 65 °C for 15 min. 

Plasma miRNA microarray hybridization 

The Hy3 labeled sample was hybridized on a miRCURYTM 
LNA Array (v. 18.0) (Exiqon, Denmark) according to array 
manual. The total 25 [tl mixture from Hy3 labeled sample was 
mixed with 25 |xl hybridization buffer, denatured at 95°C for 
2 min, incubated on ice for 2 min, and then hybridized to the 
miRCURYTM microarray for 16-20 hours at 56°C in a phalan 
hybridization bag. 

Microarray scanning and image data processing 

A GenePix 4000B microarray scanner (Axon Instruments, 
Foster City, CA) was used to scan the microarray, and collect 
fluorescence intensity. Genepix Pro 6.0 software (Axon) was used 
for image analysis. The original signal of each probe was corrected 
by subtracting that of the background. Probe intensities were 
corrected by subtracting background intensity and normalized 
using median array normalization, where each background- 
subtracted intensity was divided by the median intensity of non- 
control probes with corrected valueSSO on the array. Then, the 
fold change in miRNA expression between VSD group and 
control group was calculated for each probe (VSD group/ control 
group). The fold change of > 1.5 was used to screen differentially 
expressed miRNA. All microarray data reported in the manuscript 
have been deposited at GEO (accession number: GSE54675). 

Plasma miRNA target gene prediction 

Three databases including targetscan (http://www.microrna. 
org/microrna/home.do), mirbase (http://pictar.mdc-berlin.de), 
and Miranda (http://www.targetscan.org/vert_60/) were used 
to predict target genes for the miRNAs. 36 differentially expressed 
miRNAs were identified by microarray analysis. 15 miRNAs 
corresponding to 447 target genes were obtained from three 
databases including targetscan, mirbase, and Miranda. These 15 
miRNAs were selected for further study. 

Correlation between miRNAs and target genes 

The 447 target genes regulated by 15 miRNA were annotated 
using the Gene Ontology (GO) database (http:/ /www. 
geneontology.org/) to obtain all functional annotations of target 
genes. Fisher's exact test was then used to calculate significance 
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Table 1. Sequences of the primers 


used in the reverse transcription and RT-PCR. 




Primer name 






Primer sequence 


hsa-miR-155-5p RT 






5'GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACACCCCTA3' 


hsa-miR-155-5p: F 






5 ■ GGGGTAATGCTAATCGTGA3 ' 


hsa-miR-155-5p: R 






5'CAGTGCGTGTCGTGGAG3' 


hsa-miR-222-3p RT 






5 ■ GTCGTATCCAGTGCGTGTCGTGG AGTCGGCAATTGCACTGGATACG ACACCCAGT3 ' 


hsa-miR-222-3p: F 






5 'GGG AGCTACATCTGGCTA3 ' 


hsa-miR-222-3p: R 






5 ■TGCGTGTCGTGGAGTC3 ' 


rno-miR-379-5p RT 






5 ■ GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCCTACG3 ' 


rno-miR-379-5p: F 






5 'GGGGTGGTAG ACTATGGAA3 ' 


rno-miR-379-5p: R 






5 'CAGTGCGTGTCGTGGAGT3 ' 


hsa-miR-409-3p RT 






5 ■ GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACAGGGGT3 ' 


hsa-miR-409-3p: F 






5'GGAATGTTGCTCGGTGA3' 


hsa-miR-409-3p: R 






5 ■ CAGTGCGTGTCGTGGA3 ' 


hsa-miR-433 RT 






5 'GTCGTATCCAGTGCGTGTCGTGG AGTCGGCAATTGCACTGGATACG ACACACCG3 ' 


hsa-miR-433: F 






5 ■ CAATATCATGATGGGCTCCT3 ' 


hsa-miR-433: R 






5 'CAGTGCGTGTCGTGGAGT3 ' 


hsa-miR-498 RT 






5 ■ GTCGTATCCAGTGCGTGTCGTGG AGTCGGCAATTGCACTGGATACG ACGAAAAAC3 ' 


hsa-nniR-498: F 






5 'GGTTTCAAGCCAGGGG3 ' 


hsa-miR-498: R 






5 ■ CAGTGCGTGTCGTGGAG3 ' 


hsa-let-7e-5p RT 






5 ■ GTCGTATCCAGTGCGTGTCGTGG AGTCGGCAATTGCACTGGATACGACAACTAT3 ' 


hsa-let-7e-5p: F 






5 ■ GGGGTG AGGTAGGAGGTTGT3 ' 


hsa-let-7e-5p: R 






5 'GTGCGTGTCGTGG AGTCG3 ' 


hsa-miR-487b RT 






5'GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACAAGTGG3' 


hsa-miR-487b: F 






5 ■ GGGG AATCGTACAGGGTCAT3 ' 


hsa-miR-487b: R 






5 ■ GTGCGTGTCGTGG AGTCG3 ' 


RT: reverse-transcription; F: forward primer; R 
doi:l 0.1 371/journal.pone.Ol 0631 B.tOOl 


reverse primer. 




Table 2. Background characteristics of the three patients with VSD and the three matched VSD-free controls analyzed with the 
miRNA microarray study. 






VSD 

patient 

D28 


Matched 

Control 

CSS 


VSD Matched VSD Matched 
patient control Patient Control 
D55 CSe Dei C82 


Gender 


male 


male 


male male female female 


Age (years) 


1.2 


1.2 


2.1 2.1 1.4 1.4 


Place of birth 


Jinan, 

Shandong, 

China 


Jinan, 

Shandong, 

China 


Dezhou, Dezhou, Taian, Talan, 
Shandong, Shandong, Shandong, Shandong, 
China China China China 


Method of birth 


Natural birth 


Natural birth 


Natural birth Natural birth Natural birth Natural birth 


History of 
taking 

contraceptive 


No 


No 


Yes Yes No No 


History of 
taking Folic 
acid 


Yes 


Yes 


Yes Yes No No 


History of CHD 


No 


No 


No No No No 



VSD, ventricular septal defect; CHD, congenital heart disease; mlRNA, microRNA. 
doi:l 0.1 371/journal.pone.Ol 0631 8.t002 
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Table 3. Clinical data for the 20 patients with VSD and the 15 VSD-free controls. 



Variables 


VSD (n = 20) 


Control (n = 15) 


P-Va\ue 


Male/female 


11/9 


9/6 


0573 


Age (years) 


1.23±0.51 


1.28±0.37 


0.354 


Gestational age (weel<s) 


38.35±2.76 


39.27±0.88 


0.121 


Birth weight (kg) 


3.19±0.50 


3.12±0.50 


0.804 


History of abortion, n (%) 


9(45.00) 


1 (6.67) 


<0.001 


History of tal<ing contraceptive, n (%) 


2(10.00) 


1(6.67) 


0.497 


History of CHD, n (%) 


0(0%) 


0(0%) 


>0.05 


Maternal age (years) 


25.3±4.86 


27.13±4.61 


0.58 


History of tal<ing Folic acid, n (%) 


17(85.00) 


14(93.33) 


0.125 



VSD, ventricular septal defect; CHD, congenital heart disease. 
doi:l 0.1 371 /journal.pone.Ol 0631 8.t003 



level (P-value) of each function, which was used to screen genes 
with significant functions. A function with a P-value^O.OS was 
considered as significant. The smaller the P value was the more 
significant the function of GO was. 

Verifying differentially expressed miRNAs with RT-PCR 

The hsa-miR-93 miRNA was used as an internal reference, to 
which the results of other microRNAs were normalized because it 
was stable and expressed similarly between samples. Primer 
sequences for all microRNAs were designed based on microRNA 
sequences from the miRBase database. 100 ng of total RNA was 
used in 20 |J,L of the reverse-transcription reaction containing 
2 JtL of 10 xreverse-transcription buffer, 0.3 jiL of RNase 
inhibitor (40 U/|xL), 2 |iL dNTP(2.5 mM each), 0.2 nL of multi- 
scribe reverse-transcriptase (200 U/|j.l, EPICENTRE Biotechnol- 
ogies, USA), 0.3 |J.L of miRNA-specific primers (1 Um) and 
adding HjO upto 20 |lL. RT was carried out using a Gene Amp 
PGR System 9700 (Applied Biosystems, USA) at 16°G for 30 min 
followed by 42°G for 40 min and then inactivated at 85°G for 
5 min. After this the reaction was placed on ice ready for use. 
Meanwhile, 2 pi of RT products was added into an 8 |J,1 real-time 
PGR assay reagents containing 5 |xL of 2 x PGR Master 
Mix(Superarray, USA), 1 [iL of primers and 2 |J,L of water. 
Real-time PGR was carried out on an Applied BioSystems 
PRISM7900 thermocycler at 95°G denaturation for 10 min, 95°G 
for 10 s, 60°G for 60 s, and 95°G for 15 s. Fluorescence intensity 
was detected. The reaction was repeated for 40 cycles, then slowly 
heated from 60°G to 99°G. The melting curve was plotted and 
data was analyzed with 2"^^^^^ method. The RT primers and the 
PGR primer sets for each microRNA amplification are shown in 
Table 1. 

Statistical analysis 

AH the data were expressed as mean ± standard deviation, and 
analyzed using SPSS 16.0 software (SPSS Inc, USA). The t test 
was used for intergroup data comparison. For the microarray data, 
after normalization, the statistical significance of differentially 
expressed miRNA was analyzed by t test. For all of the analyses, a 
probability smaller than 0.05 (P<0.05) was considered statistically 
significant. P values were corrected for multiple testing by false 
discovery rate (FDR). 



Results 

The background characteristics of the VSD patients and VSD- 
free controls used in the miRNA microarray are shown in Table 2. 
The three VSD patients were matched in terms of gender, age, 
place of birth, method of birth, history of taking contraceptives, 
and history of GHD with those of VSD-free controls. The 
background characteristics of the 20 patients with VSD and the 1 5 
age- and gender-matched VSD-free controls are summarized in 
Table 3. The participants in the two groups were significantiy 
different in the history of abortion, but not in the other variables 
(P>0.05). 

Plasma miRNA expression determined by miRNA 
microarray analysis 

The miRNA expression patterns in the plasma of the three 
patients with VSD and the three VSD-free controls were detected 
using miRNA microarray analysis by miRGURY LNA Array 
(version 1 1 .0) to preliminarily identify miRNA expression differ- 
ences between the patients with VSD and the VSD-free controls. 
Results showed that compared with the VSD-free controls, the 
patients with VSD had 36 diflferentially expressed miRNAs (P< 
0.05) with 15 upregulated in their expression and 21 downregu- 
lated in their expression (Table 4, Fig. 1). 

Target genes of 15 in 36 miRNAs were predicted in three 
databases 

The target genes of the 36 differentially expressed miRNAs in 
the VSD and control groups were predicted using three databases 
including targetscan, mirbase, and Miranda. Among the 36 
miRNAs, 1 5 miRNAs corresponded with target genes across three 
databases. These 15 miRNAs were selected for further study 
(Table 4, boldfaced words). 

Expression validation by quantitative RT-PCR 

Qviantitative RT-PGR assays were performed to validate the 
different expression levels of these 15 miRNAs in plasma of 20 
VGD patients and 15 controls. The results showed that 7 of the 15 
miRNAs were not differentially expressed between the VGD 
patients and the controls; the remaining 8 miRNAs were 
differentially expressed between VGD samples and control 
samples, including seven downregulated miRNAs (hsa-let-7e-5p, 
hsa-miR-155-5p, hsa-miR-222-3p, hsa-miR-379-5p, hsa-miR- 
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hsa-miR-191-5p 



Figure 1. Heat map of microRNA (miRNA) microarray expres- 
sion data from plasma samples of patients with ventricular 
septal defect (VSD) (n = 3) and VSD-free controls (n = 3). Sample 
miRNA species are shown at the top. The miRNA clustering tree is on 
the bottom. The miRNA species are shown on the right. Cluster analysis 
classified the samples into groups based on miRNA expression levels in 
each sample. The dendrogram shows significantly different expression 
levels of mlRNAs among samples. Red indicates high expression of 
miRNA, and green indicates relatively low expression of miRNA. 
doi:1 0.1 371/journal.pone.01 0631 S.gOOl 

409-3p, hsa-miR-433, hsa-miR-487b) and one upregulated 
miRNA (hsa-miR-498)(Fig. 2). Table 5 shows representatives of 
target genes, mainly those involved in cardiac development 
(HANDl, NKX2-4, TBX-1, MAP2K4, GATA4, NOTCH 1, 
ZFPM2, and etc.), of the 8 differentially expressed miRNAs 
validated by RT-PCR. 

GO enrichment analysis for target genes 

The target genes predicted from patients with VSD may be 
involved in the pathogenesis of the VSD. GO enrichment analysis 
revealed that these differentially expressed miRNAs were enriched 
in cardiac right ventricle morphogenesis, and the top ten biological 
processes are shown in Table 6 and Fig. 3. The target genes 
NOTCH 1, HANDl, ZFPM2, and GATA3 had the most enriched 
GO functions, and were related with cardiac right ventricle 
morphogenesis. They were regulated by hsa-let-7e-5p, hsa-miR- 
222-3p and hsa-miR-433, which were downregulated in the 
plasma of patients with VSD (P<0.05) (Table 6). 

Discussion 

In order to evaluate the roles of miRNAs in the development of 
VSD this study screened differentially expressed miRNAs from 



patients with VSD using miRNA microarray analysis, preliminar- 
ily identified target genes regulated by these miRNAs, and 
provided important information for clarifying the pathogenesis 
of VSD at the molecular level. 

In this study, we detected plasma miRNAs in three patients with 
VSD and three VSD-free controls by miRNA microarray analysis 
and found 36 differentially expressed miRNAs in patients with 
VSD with 21 downregulated miRNAs and 15 upregulated 
miRNAs (Table 4). Three miRNA target gene databases including 
targetscan, mirbase, and Miranda were used to predict target 
genes for all diflTerentially expressed miRNAs. Among these 36 
miRNAs, 15 miRNAs corresponded to 447 target genes were of 
general consensus across the three databases and selected for 
further study. Expression regulations of 8 miRNAs among these 1 5 
miRNAs were validated by real time PGR: hsa-miR-409-3p, hsa- 
miR-487b, hsa-miR- 155-5p, hsa-miR-379-5p, hsa-miR-222-3p, 
hsa-let-7e-5p and hsa-miR-433 were downregulated, while has- 
miR-498 was upregulated, showing the two methods gave 
consistent results. Among target genes of these 8 miRNAs, some 
were associated with cardiac development. Further functional 
significance analysis of these target genes showed that those 
associated with the most enriched GO functions regulated 
development of right ventricular functional morphology. Key 
target genes included NOTCH 1, HANDl, ZFPM2 and GATA3, 
which were regulated by three miRNAs of the eight miRNAs 
validated by RT-PCR: hsa-let-7e-5p, hsa-miR-222-3p and hsa- 
miR-433. The expression of the three miRNAs were significantly 
(P<0.05) downregulated in the plasma of patients with VSD. 

The results presented here suggested that hsa-miR-433 may 
regulate target genes NOTCHl and GATA3. Previous studies 
have shown that NOTCHl is closely related with ventricular 
development [22]. During embryonic development, NOTCHl is 
expressed in the endocardium at the bottom of the trabecular 
layer. NOTCHl knockout in rats causes trabecular defects [21], 
and NOTCHl expression is downregulated in fetal myocardial 
tissue with VSD [23]. 

The GATA family has six genes (GATAl-6). GATA4 is a 
transcription factor that plays a regulatory role in a heart-mediated 
regulatory network, and is involved in cardiac development and 
expression of functional genes [24] . The GATA4 gene plays a key 
role in compartmentalization of the heart, atrial and ventricular 
development, atrioventricular valve formation, and arterial trunk 
separation, and it is essential for the heart to develop, mature, and 
play its physiological functions [25]. GATA4 may also be 
important for repair and reconstruction post-myocardial damage 
[26]. GATA3 is specifically expressed in T helper 2 (Th2) cells, 
and is important for Th2 cell proliferation, differentiation and 
development [27]. GATA3 gene expression is significantly lower 
in children with a combination of VSD and pulmonary arterial 
hypertension (PAH) than in those with PAH, so low GATA3 gene 
expression may be associated with VSD occurrence [28]. 

The hsa-miR-222-3p miRNA is involved in the regulation of 
the target gene ZFPM2, also known as transcription factor 
GATA4 (FOG2), another member of the GATA family. FOG2 is 
essential in heart development and can specifically interact with 
the GATA4 gene, inhibiting GATA4-dependent transcription 
[29]. FOG2 knockout mice die from heart defects on day El 3.5 of 
the embryonic period [30], hence FOG2 may be associated with 
the occurrence of VSD. 

Another of the four target genes predicted in our study, 
HANDl, is an important member of the HAND family. As a 
transcription factor regulating heart development, HAND2 is 
closely associated with ventricular muscle cell differentiation, the 
conduction bundle, the cardiac outflow tract, and epicardial and 
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Table 5. Target genes associated with cardiac development predicted for 8 differentially expressed miRNAs in patients with VSD. 



miRNA 


Genes 


hsa-let-7e-5p 
(downregulated) 


HANOI, ABCB9, ACTAl, ANKRA2, ANKRD49, BRF2, CLPl, NKDl, NME6, TRIM71, XKR8, ZNF583, NEK3, AURKB, SLC12A9, SMUGl, ADRB2, 
SLC35D2, UFMl, EZH2, TIMM17B, GPR137, ZNF354A, ATP2A2, NME4, BZW2, LIPH, DTX2, BACE2, DHX57, TARBP2, GALE, MDFI, RPUSD3, CCL7, 
0LFM4, TRAPPC1, KIAA1539, RNF20, KLHDC8B, FAM103A1, PRTG, BZWl, GNPTAB, COIL, AP1S1, HIF3A, TTLL6, TNFSF9, PLEKHG6, GNG5, 
C8orf58, GATM, B3GNT1, CTHRCl, TGDS, GATA4, RTKN, USP21, LRIG2, AVEN, PRSS22, ZC3H3, KLHL6, ABCC10, YTHDF3, STXBP5, PLEKHOl, 
nnviQR Mim ii i3 paci r" dcodci/ — j mtdi^^ ^ylPmQ di rr\r'iA ^firATAA cdatat \ r\D kid i Dif~3 arct^. dhi i NincTi mu'MiA 

L/Ua 1 yo, INIUZ, ILI,3, tAjLu, rjUnjIl-Z, IM 1 nl\5, IvItUZcS, rLUJ, \^\J\^dh, iVluA I ^A, jtA I Az, LUK, IjLK, Lnlo,3, AKl_L.D, rULL, WUdiZ, LUMMIA, 

CRCTl, CTNS, THRSP, CHRD, ZNF341, KCTD16, QARS, PQLC2, ESPLl, SLC20A1 


hsa-miR-155-5p 
(downregulated) 


ACTAl, ACTRIO, BRD1, CHD7, CSFIR, DETl, FBXOl 1, IL13, LAMP2, MYOlO, PHC2, WEEl, EOS, DHX40, HIVEP2, RNF149, TERFl, LSM14A, VPS18, 
<;^ylAD("A/1 cai i i Dmpi i nvMr"iii rrtt — 73 raiadti 1 uiqp^a dmpit^ Air~na tdv 1 i^iaai^ii; cdiw^ iddiq uiddi tdtb^q cii^iTAf; 

jIVIAnt-A't, J ALL 1 , Dl_UKL 1 , U Y InL II 1 , \^U\-/ 5, DAIArZL 1, n^r^A, KINrlzj, AH^L/A, 1 DA- 1 , MA A \ / \ D, irllNo, Lnr ID, Mori, Z.dId^o, I Zf\tD, 

USP43, STXBP5L, SAP30L, CSNK1G2, CARHSPl, TPRKB, JARID2, SGK3, PDE7A, ARRB2, C8orf4, SOCSl, USP8, KIAA1267, BOC, TRIM32, CLCN5, 
PTPN2, MYLK, MAP3K10, PCDH9, C0L21A1, PSIP1, SDCBP, TAPTl, DNAJB7, WDR45, CEBPB, MBNL3, SKIV2L2, SPII, DCLREIA, FAM105A, IKBKE, 
SHOX, RCN2, MGP, AGTRAP, GNAS 


hsa-nniR-222-3p 
(downregulated) 


ANKRDIO, CDKNIB, FRAT2, MAPKIO, NTF3, PCMTDl, RIT2, TOX, YTHDCl, ZFPM2, CYR61, FNDC3A, MIDN, KSRl, RSBNIL, SNCB, GRM1, MIA3, 
NSMCE4A, PANK3, CCDC64, ADAMTS6, GNAI2, KHDRBS2, FOS, MESDCl, RECK, VGLL4, PLCL2, KIAA1370, TBX1, HOXCIO, PAIP2, LRFN2, KRT81, 
ZFYVE16, PCDHA4, PBX3, TP53BP2, IRX5, DMRT3, HEXIMl, MBD2, TCF7L2, RBM24, SEMA3B, ZNF181, AP3B2, SMARCA5, CDKNIC 


hsa-nniR-379-5p 
(downregulated) 


CCNB1, CDKN2AIP, EIF4G2, IL28RA, KLHL14, MTMI, PCGF5, PXTl, SLC20A1, YARS, ZBTB26AKT1, GLISl, NKX2-4, PDHX, T0P2B, UTP11L, 
ZDHHC20, UPKIB, GPBPl, MSH6, CWF19L1, DENND2C, SH3BGRL3, SNX13, NDN, LM04, AP0BEC2, NR4A2, EWSRl, UBE2D2, REPSl, LHX2, ATF7, 
PAM 


hsa-miR-433 
(downregulated) 


ABCF2, CDC27, C0X6B1, GATA3, H0XA5, NOTCHl, IGFBPl, XRCC6BP1, GEMIN5, CFTR, SMG5, CLYBL, HIVEPl, RECQL, ATP5S, KIAA0195, TIPRL, 
EPHA5, C0Q3, GPATCH8, NR2F6, SLC17A5, CENPJ, CHMPS, PPMIA, HIVEP2, TDRD5, SYNCRIP, C0X8A, ATP6V0A1, NMT2, E2F3, HP1BP3, 
TBC1D19, ACAD8, RAD9A, FBX05, TRIP12, PCCB, B4GALT3, SRFBPl, ADRAIA, NME5, UBR1, LRRTM2, PAK4, STMN4, WDR45L, ZC3H6 


hsa-miR-487b 
(downregulated) 


MAP2K4, NELF, ZNF219 


hsa-miR-498 
(upregulated) 


ABTBl, ANKRDll, CRLSl, JPH3, SMAD4, TEX261, TNNTl, ZEB2, KIAA1539, MAGEEl, CRH, C5, WDR47, COLlAl, MTDH, SPAG16, TIMM17A, 
GALNT7, PCMTDl, TRIM63, ITSNl, PAM, DCBLD2, LPP, C9orf5, BECNl, RPRM, ST8SIA2, TJP2, C10orf76, JHDM1D, SACM1L, DEDD, SLC25A12, 
NFIB, MGA, CACNAIA, USP34 



miRNA, microRNA; VSD, ventricular septal defect. 
doi:l 0.1 371 /journal.pone.Ol 0631 8.t005 



endometrial development [31]. HANDl is combined with 
GATA4 in four binding sites of the promoter region during the 
regulation, both playing a regulatory role [32]. HANDl may be 
involved in the occurrence of atrioventricular septal defects [33]. 
miR-1 HAND2 inhibition results in normal development of the 
heart [34]. However, little information is found on miRNA 



regulating HANDl and VSD. In our study, we found that hsa-let- 
7e-5p was significandy downregulated in the patients with VSD. 
Compared with the control group, the hsa-let-7e-5p in the patients 
with VSD was downregulated by 0.27-fold (P = 0.024) (Fig. 3). 
This indicated that hsa-let-7e-5p may play a regulatory role in the 
development of VSD. The differentially expressed miRNAs might 
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Figure 2. Independent validation of the differential expression of microRNA (miRNA) by reverse transcriptase polymerase chain 
reaction (RT-PCR) in ventricular septal defect (VSD) patients (D, N = 20) and controls (C, N = 15). hsa let-7e-5 (P = 0.024),hsa-miR-155-5p 
(P = 0.033),hsa-miR-222-3p (P<0.001),hsa-miR-379-5p (P = 0.013),hsa-miR-409-3p (P<0.001),lisa-miR-433 (P<0.001),hsa-miR-487b (P = 0.025),hsa-miR- 
498 (P<0.001). The relative expression of eight miRNAs was normalized to expression of tlie internal control (has-miR-93), using the 2^^^*"^ method. 
doi:1 0.1 371/journal.pone.01 0631 8.g002 
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Table 6. Gene ontology analysis of the target genes of the mlRNAs. 





Term 


Count 


P-value 


Genes 


cardiac right ventricle morphogenesis 


4 


0.0002 


NOTCHl, HANOI, ZFPM2, GATA3 


respiratory system process 


4 


0.0001 


ADRB2, PBX3, H0XA5, NDN 


GTP biosynthetic process 


3 


0.0029 


NME4, NME5, NME6 


DTP biosynthetic process 


3 


0.0029 


NME4, NME5, NME6 


UTP metabolic process 


3 


0.0029 


NME4, NME5, NME6 


CTP biosynthetic process 


3 


0.0037 


NME4, NME5, NME6 


pyrimidine ribonucleoside 
triphosphate metabolic process 


3 


0.0037 


NME4, NME5, NME6 


pyrimidine ribonucleoside 
triphosphate biosynthetic process 


3 


0.0037 


NME4, NME5, NME6 


CTP metabolic process 


3 


0.0037 


NME4, NME5, NME6 


pyrimidine nucleoside triphosphate 
biosynthetic process 


3 


0.0047 


NME4, NME5, NME5 



*the table only shows the top ten gene ontology (GO) terms in the functional significance list. 
mlRNA, microRNA; GTP, guanosine triphosphate; UTP, uridine triphosphate; CTP, cytidine triphosphate. 
doi:l 0.1 371/journal.pone.010631 8.t006 



bind with specific transcription factors associated with right 
ventricular morphogenesis, resulting in tlie VSD, but this would 
need to be further identified by luciferase reporter gene assays. 

As a new research area, miRNA regulatory mechanisms need to 
be explored further. It is generally considered that miRNA 
regulates gene expression by repressing translation or directing 
sequence-specific degradation of complementary mRNA. Inter- 
estingly lower expression of NOTCHl and GATA3 is associated 
with VSD and yet we report reduced miRNA expression, which 
would in a simphstic model be expected to increase expression of 



these genes. It has also been suggested that miRNA may also 
function to induce gene expression for example, miR-373 and its 
premature hairpin (pre-miR-373), induced E-cadherin and 
CSDC2 expression by targeting specific sites in gene promoters 
[35]. In that case the target genes would not have been identified 
in this study. 

While these miRNA target genes suggest that the miRNAs may 
be involved in the development of VSD, it is worth considering 
that their altered levels may also be related to secondary 
consequences of the VSD. This is suggested because VSDs 



Sig GO terms of DE gene-BP 



cardiac right ventricle morphogenesis 
UTP metabolic process 
UTP biosynthetic process 
GTP biosynthetic process 
respiratory system process 
CTP metabolic process 
pyrimidine ribonucleoside triphosphate biosynthetic process 
pyrimidine ribonucleoside triphosphate metabolic process 
CTP biosynthetic process 
positive regulation of osteoclast differentiation 
pyrimidine nucleoside triphosphate biosynthetic process 




ID 



12 



Fold Enrichment ((Count/Pop. Hits)/(Ust.Total/Pop.Total)) 



Figure 3. The microRNA (miRNA) Gene Ontology (GO) analysis. The blue bars represent the fold enrichment value, which is a ratio of Count/ 
Pop. Hits and List.Total/Pop.Total. The P-value was calculated by Fisher's exact test. The lower the P-value is, the more significant the GO Term is (P< 
0.05 is recommended).The team of cardiac right ventricle morphogenesis's fold enrichment value is 12.23 (P = 0.000235673). 
doi:1 0.1 371 /journal.pone.01 063 1 8.g003 
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develop early in the embryo [7], so the subjects of this study with 
an average age around one year will already have developed their 
VSD. Nevertheless, whatever the reason for their altered levels in 
the VSD patients these miRNAs may prove important as 
diagnostic markers for VSD. 

Altered miRNA levels are common in a wide range of diseases 
and physiological processes from cancer to metabolism [8] . They 
have been suggested as tumor markers for gastric cancer [36,37] 
and have been implicated in the development of CHD [38-40]. 
Considering the wide range of both pathological and physiological 
processes miRNAs are involved in [8] it is possible that the 
expression profiles of the miRNAs shown in this study are the 
result of other factors within the patient, these could be secondary 
results of the VSD as mentioned above or indeed other unrelated 
processes. To rule these factors out requires further study. 

Most CHD are polygenic diseases. Many miRNA-mediated 
genes are involved in the development of CHD, and form 
regulator)' networks via their interactions [41], affecting heart 
development. This study screened differentially expressed miR- 
NAs in the plasma of patients with VSD, and firsdy reported that 
hsa-let-7e-5p, hsa-miR-222-3p and hsa-miR-433 may cause 
cardiac abnormalities by binding with heart-specific target genes. 
These target genes including NOTCH 1, HANDl, ZFPM2 and 
GATA3 may mutually interact in the development of right 
ventricle morphogenesis, resulting in VSD. We studied here the 
heart development-associated specific miRNAs by detecting the 
functional significance of target genes and using differentially 
expressed miRNA as a non-invasive diagnostic marker, and have 
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provided a new way of thinking for further study of the etiology, 
early diagnosis, and treatment of VSD. 

This investigation has some limitations. This study needs to be 
replicated in a larger, independent population to confirm the 
utility of these microRNAs as a diagnostic marker for VSD. The 
link between the differential expression of the miRNA and the 
target genes needs be investigated further, to provide definite 
evidence of their roles in development of VSD. The role of these 
miRNAs as potential diagnostic markers for VSD should be 
considered preliminary. We plan to evaluate the levels of these 
miRNAs in cardiac tissue from VSD repair, to ensure that they are 
directiy related to VSD because detection in the plasma alone 
suggests that, potentially, they could be involved in other processes 
within the body. 

Conclusion 

The circulating miRNA profiles of three VSD patients showed 
that 7 miRNAs were downregulated and 1 upregulated when 
matched to VSD-free controls. Analysis revealed target genes 
involved in cardiac development were probably regulated by these 
miRNAs. 
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